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ABSTRACT 



The problem of collisions of Shockwaves in gravity is well known and has been stud- 
ied extensively in the literature. Recently, the interest in this area has been revived 
trough the anti-de-Sitter space/ Conformal Field Theory correspondence (AdS/CFT) 
with the difference that in this case the background geometry is Anti de Sitter in five 
dimensions. In a recent project that we have completed in the context of AdS/CFT, 
we have gained insight in the problem of Shockwaves and our goal in this work is 
to apply the technique we have developed in order to take some farther steps in the 
direction of Shockwaves collisions in ordinary gravity. In the current project, each of 
the Shockwaves correspond to a point-like Stress-Energy tensor that moves with the 
speed of light while the collision is asymmetric and involves an impact parameter (b). 
Our method is to expand the metric (g^iu) in the background of flat space-time in the 
presence of the two Shockwaves and compute corrections that satisfy causal boundary 
conditions taking into account back-reactions of the Stress- Energy tensor of the two 
point-like particles. Therefore, using Einstein's equations we predict the future of 
space-time using the fact that we know the past geometry. Our solution respects 
causality as expected but this casual dependence takes place in an intuitive way. In 
particular, g^i^ at any given point r on the transverse plane at fixed r evolves accord- 
ing from whether the propagation from the center of each of the Shockwaves or from 
both Shockwaves has enough proper time (r) to reach the point under consideration 



or not. Simultaneously around the center of each shockwave, the future metric devel- 
ops a 5-function profile with radius r; therefore this profile expands outwards from 
the centers (of the Shockwaves) with the speed of light. Finally, we discuss the case 
of the zero impact parameter collision which results to the violation of conservation 
and we argue that this might be a signal for the formation of a black hole. 
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CHAPTER 1 
INTRODUCTION 

The problem of the colhsion of two shoockwaves in four dimensions (with Lorentz 
signature) that result from boosting two black holes to the speed of light is well known 
and has been extensively studied in the literature. References [I1121E1I11E1E1EIE1IS] 
provide only a subset of the work dedicated to the particular problem (for some 
interesting recent developments see [lOl [H]). In particular, [H |2l [3] deals with this 
problem pertubatively and in a series of three papers the authors compute the metric 
and derive a formula for the gravitational radiation. The collision is axisymmetric 
while the bulk matter that creates the Shockwaves as well as the back-reaction effects 
are not taken into account. 

Recently, the high energy physics community has also expressed special inter- 
est in the topic of Shockwaves collisions in the framework of General Relativity 
II21[I31[Il[I51[ISl[IIl[ISl[ISll2n] through the anti-de-Sitter space/Conformal Field 
Theory correspondence (AdS/CFT) [HI El ESI [21 ES] • The AdS/CFT duality, allows 
to formulate a process that is associated with non-abelian gauge theories at strong 
coupling as a purely gravitational problem. As a result, one may map the problem 
of heavy ion collisions (in four dimensions) onto Shockwaves collisions in gravity in 
five dimensions. Recently, we have been involved in such problems [H [131 [HI [IS] in 



gauge theories at strong coupling whose dual five dimensional gravitational descrip- 
tion, exhibits similar characteristics with those of shockwave collisions in (ordinary) 
four dimensional gravity. 

Our goal in this work is to attempt to apply our earlier experience [121 [131 [IH 
[15] and especially apply the technique that we have developed in [12] in (ordinary 
gravity in) four dimensions. We choose to work in a different coordinate system than 
[H m [3] G in order to retain the geometrical insight of the collision and in addition 
we take into account the matter responsible for the creation of the Shockwaves and 
the back-reaction effects as well. Furthermore, the collision we consider here is not 
axisymmetric but is involves a non zero impact parameter. 

We organize the paper as follows. 

In chapter [2] we state the problem we want to solve and construct the main set up. 
Our goal is to determine the evolution of the geometry assuming that we know it in 
some time interval (negative times). The initial geometry is given by two Shockwaves 
which correspond to a non zero Stress- Energy tensor Q Our method is to construct a 
perturbative approach by expanding the metric around the background given by the 
flat metric along with the two Shockwaves. Equation (12. lip shows the form of the 
metric at all times while flgure 12.31 offers a diagrammatical intuition of the terms of 
the metric we attempt to calculate in this project. 

In chapter [3] we take into account the interaction of the one particle with the 



""^More precisely our coordinate system coincides with the initial coordinate system of [1] [H |3]. 
However, in the process of the calculation [11 [21 13] choose to change coordinates. 

^This Stress-Energy tensor is due to two point-like particles moving opposite to each other (see 
figure 13. ip with the speed of light and begin to interact for positive times 



gravitational field created from the other and vice versa. In terms of Feynman dia- 
grams, loosely speaking, these corrections correspond to the diagrams of figure 13.21 
The corrections to the Stress-Energy tensor corresponding to these diagrams along 
with the corrections of the metric tensor corresponding to diagram of figure 12.31 form 
a consistent set u of the corrections that have to been taken into account. We verify 
that the modified Stress-Energy tensor is conserved (to the order of the expansion 
that we are working) and we find that it is traceless. This last condition results into 
some pleasing simplifications for Einstein's equations (compare (12. 5p with (I3.22p ). 

Chapter H] deals with the field equations and the specification of the gauge. Our 
attempt is to perform a perturbative calculation about a metric that looks almost 
fiat but contains two Shockwaves moving opposite to each other and colliding (see 
(12.111) ). The Shockwaves provide an effective Stress-Energy tensor in addition to 
the (actual) Stress- Energy tensor (see (14. 3p ) that creates the two Shockwaves; these 
terms correspond to products of the form tit2 and T^,y of equation (14. 4p respectively. 
A suitable gauge choice simplifies the field equations (14. 4 p to equations (14. 9 p which 
are solved in the next chapter. 

In Chapter |5] we specify the boundary conditions of the field equations (14. 9 p and 
the corresponding Green's function. In particular, we seek for causal solutions and 
therefore the associated Green's function to the differential operator (I4.10p is the 
retarded Green's function. We show how the integrations on the light-cone and 
transverse plane may be performed omitting some of the intermediate steps for Ap- 
pendices 13 E] and O Eventually we derive a formula for the metric tensor, equation 
(IS.lOp . to the order we are working. This is our final result and it generally has the 
structure (16. 2p . 

Finally in Chapter [6] we discuss about the area of validity of our solution and 



By consistent we mean with the order of the parameter we are expanding (see also (|4.2|) 
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summarize our conclusions. In particular, we argue that the presence of matter 
and the back-reaction effects may not be ignored as they result to an important 
contribution to the metric. We also see that as the impact parameter tends to zero, 
the metric diverges logarithmically. It is in our belief that this is a signal that a 
classical approach to the problem stops being valid and that a quantum description 
is required. Lastly, we talk about the general form of the metric. Although its 
evolution is constrained by casualty as expected; this evolution takes place in an 
intuitive way: At a given (proper) time, any arbitrary point on the transverse plane 
evolves according from whether the signal from the center of the one or the other 
Shockwave or both, has enough time to reach the point under consideration. We had 
encountered such a behavior in an analogous set up in [12] although the geometry 
there was Anti de Sitter geometry in five dimensions. In [T2] we had claimed that 
a similar evolution of the metric should also be observed in four dimensions; in this 
work we verify our conjecture. 



CHAPTER 2 
SETTING UP THE PROBLEM 



2.1 Single Shockwave Solution 

We begin by defining the coordinate system (gauge) we work. We choose to work in 
hght-cone coordinates defined by 

x^ = {x~^ , x~ , x^ , x'^) x^ = -^— (2.1) 

V2 

where x° is the time axis and x^,x'^,x'^ cover M'^. The convection for the fiat metric 
that we use is 

g^, = {-1,1,1,1) (2.2) 

We suppose that we have a black hole metric that we boost to the speed of light 
along a given direction (x^ direction). The metric is known |4] and is given by 

ds'^ = g^i^dx^dx'^ = —2dx^ dx~ + ti{x^,x^,x'^)dx~^'^ + dx^, (2.3) 

where 



t^ = -i_iS{x^)\og{kr) r = (x\x^) r = ^/{x^y + (x^)^. (2.4) 

According to equation (12. 3p we denote the transverse fiat metric by dx^ = (dx^Y + 
[dx^Y while 5{x^) denotes the delta (Dirac) function. The parameter jj, has di- 
mensions of length and its physical meaning will become apparent in what follows 

5 



while k serves as an ultraviolet cutoff and whose physical meaning is discussed in the 
conclusions (see Section ( 16.21) '). 

n 

One may check directly whether (12. 3p solves Einstein's equations u which may be 

CclSTJ cLS 

R,u = ^l (t,, -^g^uT] T = T; = T^, ^^^ 4 = 8776-4 (2.5) 

where R'^'^ is the Ricci tensor, T^^, the Stress-Energy tensor and G^ the Newton's 
constant (in four dimensions). 

Direct substitution of (12. 3 p in the formula that computes R^^, results to 

i?,. = V^.+ (y/i5(x+)5(2)(f)) (2.6) 

where 5^+ is a Kronecker delta. This implies that all components of R^^ are zero 
except from -R++. In order to arrive to equation (12. 6p we had to evaluate the following 
linear differential expression 

R++ = -lvlt,{x^,x\x'). (2.7) 

where V^ is the Laplace operator in two dimensions while we used the identity 

Vl\og{kr)=27r6'-^\f). (2.8) 

Equations (12.51) and (12. 6p imply that the metric tensor of equation (12. Sp corre- 
sponds to a Stress- Energy tensor. Shifting the origin □ along negative x^ for distance 
b we find that the Stress-Energy tensor is given by 



""^The cosmological constant is assumed zero. 
The reason for this shift will become apparent in the next section. 



' /ilogC'T) 



Figure 2.1: The shockwave solution: a graviton (curly line) is emitted from the bulk 
source (straight line) with coupling nlog{kr) which is measured at point 
x^. This is a very special case where a single graviton exchange between 
the source and the bulk happens to be an exact solution to the non linear 
Einstein's equations. 



T« = V5.+ (^6{x^)6{x' - b)6{x')] . (2.9) 



The presence of the superscript ('^-') on T++ of (12. 9p is to highlight that it is of first 
order in the parameter /i q This equation implies that the shockwave of (12. 3p is a 
consequence of a point particle moving along the x~ direction with the speed of light 
and hence its massless. Indeed, the ratio — has dimensions of mass as should. One 
may check that T^^, of (12. 9p is covariantly conserved. 

We may gain some insight in the physical system at hand if we use a diagram- 
matical approach. Despite that General Relativity is a non linear theory, the metric 
(12. 3p satisfies a linear differential equation. The fact that (12. 3 p is an expression of 
a perturbation of a flat metric proportional to fi suggests the diagram of figure 12.11 
It represents the measurement of the gravitational field at point x^, which, loosely 
speaking, is created by a single graviton emission from the source (point-like Stress 



Generally, the superscripts emphasize the number of times the source ti (i=l,2 assuming we have 
two sources) appears. Generally in the object A^"' there exists the product ^1^2^ or any linear 
combination of differentiation/integration of ii and ^2 with respect to their arguments x^. 



Energy tensor) of equation (12. 9p with effective coupling proportional to fi\og{kr). In 
other words for this special case of a single shockwave, the first order solution happens 
to be the exact solution to all orders. 

2.2 Superimposing two Shockwaves 

Having defined all the necessary ingredients we now proceed to the main part of the 
setup. We want to superimpose two such Shockwaves whose sources are two point-like 
distributions of matter moving towards each other U in space-time. We want to collide 
these Shockwaves (and as a result the corresponding Stress-Energy tensors as well) at 
a non-zero impact parameter and hence study the problem within the classical theory 
of gravity. Therefore, T^^, has in addition to (12. 9 p the symmetric part 

yW ^ ^s{x-) S{x' + b)S{x^) (2.10) 

which creates a second shockwave. In terms of space-time, this would correspond in 
"colliding" two metrics in an off center process. Figure 12^2] represents the four dimen- 
sional picture, right before the collision of the two Shockwaves. Following |13j.|15]. 
the metric that describes the process should look like 

ds'^ = —2dx~^ dx~ + dx\ + t\ [x'^ ^x^ — 6, x^) dx'^'^ 

+ 4^^(x-,x^ + 6,x2) dx-^ + 9{x+)9{x-)gj^J{x^ , z)dx''dx'' + ..., 
t^ll{x^Tb,x^) = -/ilog (k^{x^ T by + {x^)A 6{x^). (2.11) 

The first three terms correspond to the flat (Minkowski) space. The next two are of 
first order in fi and are created by the two point-like particles. These move (initially) 

This is true for negative times only. For positive times their trajectories are ahered and have to 
be specified taking into account back- reaction effects (see Chapter [3]). 
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fi6{x } log■(^■r2 




fi6{x+)log{kri) 



Figure 2.2: The two Shockwaves before and after the colhsion moving along x axis 
and dragging a perpendicular gravitational field which is constant along 
the circular lines. They collide at the origin producing a gravitational 
field in the forward light cone. Our goal to compute the "produced" 



metric and in particular g 



(2) 



' ii\og{kri) 




■ fJ.log(krn) 



Figure 2.3: The Feynman diagram that represents the /i^ correction of the metric: 
It represents along with the diagram of figure 13.21 (see Chapter [3l), the 
first non trivial correction to (12.111) . It shows how the two metrics that 
each one looks like (12. 4p merge. The gravitational field is measured at 
the point x^. 



towards each other along x^ and they have an impact parameter 26 along the x^ 
axis as figure 13.11 depicts. As they are the sources of the two Shockwaves, they 
correspond to two vertex diagrams that look like the one in figure 12.11 This is a 
superposition of two metrics with each one looking like (12. 4p . However, the non- 
linearities of the gravitational field require higher order terms. The second order 
corrections are explicitly displayed in (12. lip and they appear once the two Shockwaves 
cross each other; in the forward light cone. This is precisely the meaning of the 9- 
functions; they emphasize that the metric (12. lip solves Einstein's equations exactly 
in the presence of both Shockwaves only for negative x^. The additional terms of the 
metric appear in the forward light cone only and describe the effects of collision. The 
main work of this paper is to show how these terms may be calculated to order /x^, 
that is find g^^f}. The second order correction in /x of g^j^^ corresponds to the diagram 
of figure 12.31 
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CHAPTER 3 
BACK-REACTIONS 

3.1 Corrections to T^^, and Geodesies 

As has been already mentioned below (12. 9p . T|^ is conserved in the gravitational field 
o f J2.3p . (12. 4p . In fact, conservation for this case happens to be valid to all orders in 



V*V« = 0. (3.1) 

Conservation to first order is still valid when we consider simultaneously T)_^ and 
T_l in the presence of the gravitational field (12. lip . However, this is no longer true 

at the second order in /i. The reason is because the T ( T++) source moves in the 

gravitational field of the ti (^2) shockwave, altering its initial trajectory. Figure 13.11 
outlines what happens while figure 13.21 offers a diagramatical intuition regarding the 
self-corrections to T^^. This implies that we should correct T^J,-* in order to preserve 
conservation (of the total Stress- Energy tensor). However, since we do not know the 
nature (equation of state) of T^^ we make the assumption that these objects interact 
only via gravitational forces Q 



In practice, only the first order in ^ appears in the resulting equations. This is in accordance 
with our intuitive picture of figure [^TTl Gravity behaves linearly with respect to the metric (|2.3p . 

More precisely, we assume that any other interactions are small compared to the gravitational 
forces. 
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t> 




t<0 



t>Q 



Figure 3.1: The sources (point-like particles) represented as black dots and their tra- 
jectories: For negative times they move along straight lines. At t=0 each 
particle (dot) intersects the shock due to the other particle and its tra- 
jectory suffers a kick (see curved path); hence T^jy changes with time. 



T_ 



++ 



n 



^^Maaap 



X'' 



^^^^j^jxAA A/^ 



r_ 



r_ 



Figure 3.2: Graviton emission diagrams resulting from self corrections to T^^,. The 
first source interacts via the gravitational field created by the other and 
vice versa. The point x^ is the space-time point where g^j^y is measured. 
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3.2 Calculating the Corrections for T^i, 

Since these particles are point-like and massless, they should travel along null geodesies; 
as it is rigorously shown in [26] conservation of (the total) T^y is then guaranteed. 
This suggests that we need 

^'^ = ^ E 4)4)7=^^'^ (%) - %)(^«)) (3.2) 

71 11 

which gives the total Stress -Energy tensor of both point-like particles of mass -^ 
each moving along the trajectory X(/)(s(/)) parameterized by S(/). The quantity g is 
the determinant of the (total) metric tensor, the factor vr is in agreement with our 
convention (that reproduces (12.91) - see below) while the dots denote differentiation 
with respect to the parameter si. Before calculating higher order corrections to T^^^, 
we find it instructive to check whether this formula reproduces (12. 9p in the case of 
one particle (/ = 1). The trajectory of this particle which moves ultra-relativistic 
ally along negative x^ is parameterized by 



a;f^)(a;°,x\a;^a;=^) = (t,0,0,-t). (3.3) 

In light cone coordinates and choosing to parameterize the trajectory by x~ , equation 
(13. 3p then implies 

xj^) = {x^{x-) = 0,x'{x'y,x\x') = 0,x^(x~) = 0) = (0,x";0,0). (3.4) 

Direct substitution of (13. 4p to (13. 2 p yields 



T-- = T(i)++ = -^S{x^)5ix')Six') (3.5) 



Ka 
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which is exactly equation fl2.9p (for 6 = 0). This computation also clarifies the 
convection of vr in (13. 2p . The part of T^,^ of equation (I2.10p due to the second particle 
is reproduced similarly. 

The next step is to calculate the next (second) order corrections (in /i) of T^^, which 
translates into finding the corrections to the trajectories X(7). These may be obtained 
from the geodesic equations. In particular, we only need the first order corrections to 
a;(7) as we already have a power of /i in front of the summation operator (see (13. 2p ). 
The geodesic equations we need read 

^'^I + ^J;up^"lX''l=^ ^,^=1,2 I^J (3.6) 

and are interpreted as the motion of particle / in the gravitational field of the particle 
J (due to Fj. where F are the Christoffel symbols) and vice versa; this is precisely 
the meaning of the subscripts / and J. We begin with computing the corrections to 
the particle 1 = 1 whose (first order) perturbed trajectory looks like u 

(x^(a;-))(°) + (x^(x-))« = ((x+)«,a;- + (a;-)^, 6 + (xi)«, (a;^)^) (3.7) 

where we have chosen to parameterize the trajectory with x~ (i.e. S(/=i) = x~). The 
superscript ^'^' denotes the order in the expansion. Taking into account (13. 2p and the 
fact that F/'2)^„ ~ /i otherwise is zero, we deduce that both of the terms x'^ or x^ have 
to be of zeroth order; i.e. the only choice is u = p = —. This implies that we need 
to determine F^_ to first order in /x that arise from the second particle u and which 
read 



F = —-t2,x- F — F — — -^2,^1 F — —-^2,1:2. (3. 



^We drop the subscript (i) which labels the particle 1=1 for simplicity. 

^Where we (also) dropped the subscript (2) which labels the second particle. 
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A few explanations about our notation are in order: the term ^2 is due to the second 
particle and is given in equation (12. lip . The subscript 2,x^^ on ^2 denotes ordinary 
differentiation of the source u ^2 with respect to the coordinate x^ . According to 
(12. lip , the source ^2 is of first order in /x |j and as a result the same applies for the 
r's of (13. 8p : it should be by now obvious that the Christoffel symbols are due to the 
second particle (1(2 )• 

The final step is to integrate (13. 6 p using (13. 8p and using causal boundary conditions 
and substitute the result in formula (13. 2p . As we are interested in second order 
corrections in /i, we immediately conclude that at least one of x^^ or x^ should be of 
order zero, i.e. /i = — or z/ 7^ — or yU 7^ — , z/ = — . We also note that ^/^ ~ H-O(yLi^) 
and hence according to (13. 2p corrections from \/^ do not contribute to T^^ at 0{fi^). 
We consider two cases 

Case I: y/^ = 1, /i = z/ = — 



In this case the modification of T^i, of equation (13. 2 p to this order we are working 
is in the arguments of the delta's. We have l 

X 5 ix"^ / dx~ / (ix~t2,x2 1 (3.9) 



As matter implies curvature and curvature (a non-flat metric) implies matter, we will often use 
these notions interchangeably while their distinction should be evident from the context. 



,2- 



From now on we drop the superscript ^ ' which denotes the order in /i from t\ ^ 

"^ All the integrations with respect to dx that follow from now on (see also Appendices 
1X1 and |B| will imply the obvious: / dx^ stands for / dx'~ , / dx~ / dx^ stands for 

dx' I dx" etc. 
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Expanding the delta's to first order in the sources we obtain 



2 
^4 



6{x+)5{x^ - b)5{x^) - - I dx-{ t26'{x+)5{x^ - b)5{x^) 



+ f dx- (t2,^iS{x+)S\x' - b)S{x^) + t2,^2Six+)S{x' - b)S'{x^)) j . (3.10) 

We may cast (the first order correction terms of the) last equation in a compact form 
by expressing it in terms of ti and t2- Using the identity 

1 



fi6{x^)6^^\r-b,,2) 



2tt 



Vlti,2 6i,2 = (±6,0) 



(3.11) 



(see ([21]) and (1211]) ) the 0{fi^) terms of fl3:T0]) take the form 



(T(i)++)(^) = {T-,-r 



AkI 



X dx { t2V5_ti x+ + V^ti x'l / dx t2 a;i + V^ti a.2 / dx t 



(-2,1- 



(3.12) 



where we restored the subscript (i) and the superscript ^^^ in order to highlight that 
this is the second order correction to T^^ of the first particle. 

Case II: \/—g = 1, /i = — , z/ 7^ — 



In this case the modification of T^y of equation (13. 2 p to this order we are working 
is in the factor x'^x~ = x'^. Combining (13. 6 p and (13. 8p one may compute x". Plugging 
this result to (13. 2 p and employing the identity (13. lip in order to write the transverse 
delta's in terms of ti yields to 



(T^^^r = - 


-4.1^^^^^^ 


(3.13a) 


iT,i)^ir = 


-{T-'f^ = -^V% j dx-t2,x. 


(3.13b) 


(T(l)+2)(^^ = 


{T-'f^ - ^\vlh 1 dx~t2,x.. 
16 


(3.13c) 



Second order corrections to the total T^^ 

The second order corrections to the stress energy tensor {T^i^^yp"^' of the first 
particle is given by equations f l3.12p and fl3.13p . The corrections {T{2)^Jiv) of the 
second particle may be found analogously and therefore the second order corrections 
to the total stress energy tensor read 

(r+_)(2) = -(r+-)(2) = i J^ (t^Vlti + hVlh) (3.14a) 

(r++)(^) = -— / dx' U2Vlti,^+ +V\ti^^i / dx~t2^^i + Vlti,^2 / dx^ta,^^) 

(3.14b) 

(r__)(2) = 1-1 1 dx+ UiVlt2,x- + vit2,x2 /"rfx+ti,,2 + vlt2,xi /"c?x+ti,,i j 

(3.14c) 

(r+i)(2) = l^Vlti /rfx-t2,.i (T_i)(2) = 7Avlt2 /f/x+ti,,i (3.14d) 



^ Ivlti fdx-t2,.2 (T_2)^'^ = 1^ 

^i J 4/1:4 

(rn)(') = (T22)(') = (Ti2)(')=0 (3.14f) 



(T+2)^'^ = T-Vlti / dx-t2,.2 (T_2)^'^ = T-Vlt2 / rfa;+ti,,2 (3.14e) 

4: Ka I 4 /€/ 



The first equality in equation (I3.14ap is not completely obvious and so we prove it 
below by considering 

= ^f /^i(T,,)« + giW'iiT-Y'' + 9fj'iiT--r 



7^(t2Vlti+t2Vlti) (3.15) 



4^4 
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where in the fourth equahty we used the fact that the Stress-Energy tensor of the 
point particles (see f l2.9p . f l2.10p and (13.1 ip ) to first order in /x may take the form 

tS = -ivlt, T« = -^Vit,. (3.16) 

Despite working mostly with (13.141) which is a compact expression, for concreteness, 
we write T^,y of (I3.14p in terms of the coordinates in order to clarify its form. Defining 

7\=f—bi f2=r — h2- (3-17) 

and employing (13. lip in (I3.14p we obtain 



2 

(T+_)(2) = '^\og{2k\h\)5{x+)8{x-) ((5(2)(r-l) + ,5(2)(r^2))) (3.18a) 



{T,,r = f|^(-' 



\og{2k\b\)5\x^)5^'^\r-{) + -^-5{x+)5\x^ - b)5{x 



+ 4/+(l2)2 ^(^")^(^^ - ^)^'(-') 



x^ + h 

(3.18b) 



(r+2)(') = (3.18d) 

The asymmetry between TJ./ and TJig is due to the fact that the impact parameter 
b has only x^ component (see (12.110 ). The rest non zero components that complete 

(I3.18P may be obtained using the discrete symmetries of the problem: T_l and T_-l 

(2) (2) 

may be obtained from T|_j_ and T|{ respectively by interchanging + o — and b o —b. 

We want to justify our claim that the TJ^^ corrections to T^^, correspond to null 

geodesies. For this we consider the line element ds"^ = ~2dx^dx'' + t2{dx~'Y + dx±^ 

of the second shockwave t2 (the analogue of (12. 3p ). For time- like distances and for 

fixed transverse position we have ds^ = = —2dx^dx^ + t2{dx^Y and integrating 



over the discontinuity due to the trajectory of the first particle (ti) along the second 
Shockwave (t2) we deduce that 

(Ax+)(^) = dx = - dx-t2 (3.19) 

The superscript ^'^' on this equation highlights the fact that the discontinuity along 
the Shockwave ti is of first order in fi and it implies that the trajectory of the second 
particle will be modified by along x~^ by (Ax^)^^'. But according to the argument 
of the first ^-function of (13.91) . equation (13.19^ is exactly equal to the shift along x~^ 
that we have already encountered from the geodesic analysis (for the first particle). 
This completes our argument. 

3.3 Conservation, Tracelessness and Field Equations 

The second order corrections to (the total) T^^ have already been calculated in the 
previous section. One, may check by a direct computation using ( IB. 21) that this T^^ 
is covariantly conserved. Explicitly this means that 

((V/^)(o) + (WY^^) ((T^.)(i) + (T,,)(2)) = 6^,VlhVlh + 0{fi') 

= 0(/x^) when 6i ^ h otherwise ~ 5±J^^\0) (3.20) 

where V denotes a covariant derivative while the superscripts denote the order in n 
while we have used the identity flB.2p . Therefore we conclude that T^^^ is conserved 
if and only if the impact parameter is not zero q This is one of our main conclusions 
in this paper. 



The zero impact parameter b causes problem in the metric as weU. As we wiU see, as 6 — >■ 0, the 
metric tensor diverges logarithmically (see section [6.21) . 
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It is also useful to compute the trace of T^,^ as it enters the field equations (see 
(12. 5p ). A short computation yields to 

T = g^'^T^, = {gn^'\T,,Y'^ + {gn^'\T,,Y'^ = + 0{fi') (3.21) 

which shows that the stress- Energy tensor is traceless to order fx^. Tracelessness is 
very convenient as it simplifies Einstein's equations which become 

R^.u = f^lT^u + 0(/i=^) kI = SnGi. (3.22) 
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CHAPTER 4 
FIELD EQUATIONS 



4.1 Field Equations to 0{f?) 

In this section we wish to write an exphcit form of fl3.22p up to order 0{jj,^). In order 
to determine these (differential) equations we take into account that the zeroth order 
terms satisfy fl3.22p trivially as R^^^^ = T^^J = while i?|^^j (resulting from the first 
order terms of (I2.1ip ) is compensated by T'^J of equations (I2.9p and O2.10p . Thus, we 
only need 

41^ = <TfJ (4.1) 

where TJ^J has already been calculated in the previous Chapter and is given by (I3.14p . 
It is crucial to state that KJ^l receives two different type of contributions: (a) The 
contribution due to the (pre)existing Shockwaves (that is due to the ti and ^2 terms 
of (12. lip ): we denote this contribution by {R]^l)ti2- (b) The contribution due to the 
(second order) corrections (in /i) of the metric (that is due to g'j^})] we denote this 
contribution by {R]^l)g- Recalling equation (12. lip that gives the form of the metric 
at all times and expanding (14. ip to 0{^'^), we expect that it should have the form 
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where {Rj^^)ti2 (and T^^') is known while {Ri^})g is what we will use in order to 
determine gj^j) (see fl4.4p ). Equation (14.21) may also be cast in the form 



{^Hu)g — '^4^/iiy i^tiiy) 



^fll/ )tl2 



(4.3) 



and view {Rj^},)ti2 ^s an effective (contribution to the tottal) Stress-Energy tensor (see 
(14. 4p and (14.91) ). Dropping, the superscripts ^'^' from g^j^} and T^^J and the superscripts 
^^' from t\ 2 for simplicity u, we find that the components of (13.221) to second order in 
H read 



:++) 



(+-) 



— 9++,x^x^ — 9++,x^x'^ + '^9+l,x+x^ + '^9+2,x+x'^ 
~ 9u,x+x- — 922,x+x- 



l^/i,T++i 



(4.4a) 



1 r 



-1) 



~ '^9+-,x'2x'i — 2(7+-,xixi + '^9+2,x-x'i + '^9+i,x-x^ — '^9++ 

+ 2fi'_2,x+a;2 + '^9-l,x+x^ " '^9ll,x+x- " '2922,x+x- + %+-, 

~ 2(7 ,x+x+ — 2tl,xi^2,xi ~ 2ti^2.2t2,x2 + tl,x+t2,x 

~ 25'+l,x2a;2 + '^9+2,x+x- ~ '^9++,x-x^ + '^9l2,x+x^ ~" '^922,x+x^ 



kIt+., 



(4.4b) 



+ '^9+-,x+x^ + 25'+i,2:+a;- — '^9-l,x+x+ + ^l,x+^2,xi 



K4T+1 



4-^+1' 



(4.4c) 



11) - 
^ ' 2 



- 9u,x2x^ + '^9l2,x^x^ - 922,x^x^ + '^9+-,x^x^ " '^9+l,x-x^ " '^9~l,x+x^ 



+ "^gil^x+x- + ^l,xi^2,a'l + tl,x^x^t2 + ^1^2,zixi 



^4^11 = 0, 



(4.4d) 



(22) - — gii^x-^x^ + 25'i2,xia;2 — 922,x^x^ + '^922,x+x' — '^9+2,x-x^ " 25f_2,3;+a:l 



12) - 

^ ^ 4L 



2fi'+-,x2x2 + ti_3;2t2,a:2 + ti^x'^x^h + tit2,x2a; 

45'+-,xix2 — 2(7+1^2,-^2 — 25f_,_2,x-xi ~ 2(7_i ,j.+, 

^2,x2^1,xi + ^l,x'2i2,a;i + 2t2il,a'ix2 + 2tit2,xia; 



A:4'T22 = 0, (4.4e) 

'^9-2,x+x^ + %12,x+x- 

(4.4f) 



^4^12 = 



We restore the superscripts that denote the order where is necessary. 
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where ti_2 are given by (12. lip and correspond to the geometry for negative times 
while the components of T^^, are given by fl3.14p . Indeed (I4.4p has the expected form 
of equation fl4.2p . The above set of the field equations has been written without 
specifying the gauge. In the next section, we will see how these may be simplified by 
making a convenient gauge choice. 

4.2 Choosing the Gauge 

We follow the standard procedure in order to solve (14.41) : We define a new coordinate 
system x'^ with respect to the old one x'^ (see (12. ip ) by 

;J'^ = a;^ + (^'^)(2) (4.5) 

where (^'^) is an arbitrary function of the (old coordinates) x'^ and is second order 
in yU. Obviously, this transformation induces a second order change in /i to gj^J but 
does not alter gj^^J and gj^J. More precisely the second order terms of the metric 
transform to 

^2^ = ^£) + e + ^S (4-6) 



where the semicolon denotes a covariant derivativ^ij. What is remarkable is that 
the field equations remain invariant under this transformation as it may be shown 
on general grounds [H [26l ETJ |28]. For concreteness, we exhibit it here for the (12) 
component, equation ( ]4.4fp . Taking into account that at the order we are working, the 
covariant derivative may be replaced by ordinary differentiation, plugging the tensor 



^In fact gj^J — gj^J — ^f_,-v + ^u■^i is exactly equal to the action of the Lie derivative acting on g^^} 
along the vector field T, i.e. C^g^i^} = ^^2) + e^^) . 
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v!^ + ^i^Jj^ into the differential part of (14.4fp and dropping (again) the superscript ^'^' 
from the ^'s for simphcity, we obtain 

- 2 (^_,:ci + ^l,x-)^^+^2 - 2 (^-,^2 + ^2,x-)^^+^i + 4 (^l,:c2 + ^2,xi),s+x- 
1 



{(4^-1 



2eH 



2^+,x2,X-,xl) + (+ ^ -)} 



+ {(4^1,x2,a;+,x- — '^^l,x+ ,x- ,x2 — 2i^l 



a;~,xT X- 



+ (1^2)} 







(4.7) 



where we used the fact that partial derivatives commute. So far, the vector field ^ 
has been arbitrary while the result of the transformation (14. 5 p on (14. 4p is just the 
relabeling g^^, — )• g^^. A convenient choice of ^i, is the one that satisfies de Bonder 
gauge 



9f,u,^--V,.u9^^ = 



(4.8) 



where t^^j, = gj^J is the flat metric. Applying this gauge to the field equations that we 
are interested, equations (14. 4p . and dropping the tilde symbol from g^i, for simplicity, 
the field equations simplify to 

1 



(++) Og, 



dx ^2 Viti x+ 



+ Vltlj.! / dx t2,a;l + V^tl ..2 / dx t2,x2 



,x^l 



(4.9a) 



(+1) 

(11) 
(12) 



1 



^9t.. = -2 (^2Vlti + taVlti) - ti^,it2,xi 



Og, 



1 

~ ^l,x2^2,a,-2 + -tl,x+t2,x- : 
1 2 /" , - 1 



1 1 

Ogi2 — -t2,x2tl,xl + ^^l,a;2i2,xi + hh/x^x'^ + ^1^2,^1x2 — 



(4.9b) 
(4.9c) 

(4.9d) 
(4.9e) 
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where we have used fl3.14p while D is the scalar operator in flat space; that is 

n = r^d^d, = -20^+8,- + Vi. (4.10) 

In the next chapter we will see how equations (14.9P may be solved imposing appro- 
priate boundary conditions. 
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CHAPTER 5 
SOLVING THE FIELD EQUATIONS AND CAUSALITY 

5.1 Green's Function and Boundary Conditions 

In this section we will show how to solve (14. 9p by seeking for causal solutions. The 
casual boundary conditions imply that the flat metric in the presence of the shock- 
waves, as can be checked, is an exact solution to Einsteins equations with a right 
hand side given by fl2.9p and fl2.10p . for negative times only. For positive times, that 
is for x^ > and x~ > 0, the second order corrections (in fi) of the metric are 
switched on as the point x'*' = x~ = is the collision point. Simultaneously, the 
initial Stress-Energy tensor of the (massless) particles that induces the shock-waves 
suffers a change (see figure [3^ and (I3.14p ) that also has to be taken into account. 

The retarded Green's function we need corresponding to the differential operator 
(14.1 op is known and in light-cone coordinates is given by 

(5.1) 

where according to (12.40 . r = (x\ x^), 9 denotes a theta (step) function while equation 
(15. ip (the retarded Green's function) satisfies 
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nG(a;^ - x'^) = S{x+ - x'+)S{x- - x'-)S'^'^\f- r'). (5.2) 

5.2 Integration over the Light-Cone Plane 

The procedure we have to follow is standard: we convolute the right hand sides of fl4.9p 
with (15.11) and integrate in all over space-time. We find it convenient to introduce 
the following notation 

G®f{x'^)= I dx'+ I dx'~ I (fr'Gixf'-x'^')f{x'^) (5.3) 



where f{x^) is any arbitrary function of x'^ while the last integral denotes integration 
in the transverse plane. We wish to perform the x integrations for all the possible 
cases that we will encounter while specifying g^^J, from (14. 9p . We organize these 
integrations in five cases while we leave the details of the calculation for Appendix [Al 

Remark 1: In addition to the five cases of the x integrations which we solve in 

Appendix |X] there exists another case that involves the evaluation of 

G«/<;.-/d.-(*,,Vi*,,+*,.,Vi<,.,) (5.4) 

that arises from (I4.9ap . Due to the complication of the calculation we evaluate this 
term in a separate Appendix ( see Appendix [HI equation (IB.3P ). 

Remark 2: The results of all of the integrations in the light-cone plane (performed 

in Appendices 1X1 and iBll are proportional to the product 9{x'^)9{x~). This implies 
that the second order corrections to g^^p appear in the forward light-cone (see figure 
lA.ip which is what we have initially demanded by seeking for a causal solution. 
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Remark 3: The right hand sides of (14 .Qp contains expressions of the form tit2 dif- 
ferentiated with respect to x^ in some fashion. According to (12. lip , these expressions 
are proportional to 6{x^)6{x^) or their derivatives. Our previous analysis has already 
taken care of the x^ integrations and so from now on, by ti^2 we will mean just the 



transverse part of ti^2: — /ilog(A/(x^ ± 6)^ + (x^)^). 

5.3 Integration over the Transverse Plane 

Having performed the x^ integrations we move to the integration over the transverse 
plane. The quantities we have to integrate have the stTnctnTeudrci^{tit2) or 9^- j(tit2) 
where a,c;i,j = 1,2. The subscript a (c) and the superscript i (j) denotes differenti- 
ation of the source t^ (tc) with respect to the space-time coordinate x* (x-*). We may 
reduce the number of the different integrals that we have to perform by working as 
follows. We firstly introduce the vectors 

bl = ibu,bi2) &2 = (&21,&22) (5.5) 

and generalize the form of (the transverse part of) ti^2 given by (12. lip to 

— * — * 

ti(r - bi) = -filog{kri) t2{r - 62) = -Ailog(A;r2) (5.6) 

where (ri^2 where defined by (13.170 . The next step is to exchange the derivatives 
acting on ti 2, that is dr^i with differentiations with respect to 6's of (15. 5p . that is with 
~f^6„,Q Finally, at the end of our calculations we take the limits 

hi^{h,Q) b2^i-b,0). (5.7) 



There is another case where we have to mtegrate terms of the form Vj^ii.2. However, according 
to p.lip these (transverse) integrations are trivial as they involve delta functions. 



^So for instance ti^x^t2^x^ takes the form 9^^jb2j(iii2)- 
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Looking equations (J4.9bp - fl4.9el) we see that they involve the product ^1^2 differenti- 
ated with respect to the transverse coordinates u. Exchanging the transverse differ- 
entiations, according to our earher discussion in this section, with derivatives with 
respect to the components of 61^2 and taking into account the transverse part of the 
Green's function, (15.11) . we see at once that we have to calculate the following integral 

-| /»oo /'27r 

J{'ri,r2,T) = / / r'dr'd(l)'5{T -r')\og{k\r' + f{\)\og{k\f' + f2\) (5.8) 

27rr Jq Jq 

where we have introduced the convenient factor 



27rr 

Now, the non trivial integration is the angular integration as the radial one be- 
comes trivial due to the ^-function. Both of the integrations are performed in Ap- 
pendix [C] and the final result reads 

J{ri,r2, t) = 6{ri - r)6'(r2 - r) Ji(ri, r2, r) + 6{t - r2)6{ri - T)j2{ri, r2, r) 

+ e{T - n)e{r2 - T)Js{n,r2, r) + e{r - r^)e{r - n)Mn,r2, r) (5.9) 

where the jT's may be found with the help of table IC.l] and equation ( 1C.4p . Equation 
(15. 9p is the last ingredient that allows us to obtain the desired solutions for equations 
(14. 9p . We display the results in the next section. 

5.4 The Formula for gj^J 

Having performed all of the integrations arising from the convolution of the right 
hand sides of (14. 9 p with the Green's function ( 15. ip . we are in a position to derive the 
final formulas for gj^J Q The ingredients that we need, have been obtained or defined 



^Where iii2 ~ log(fcri) log(fcr2); see (|2.11|) while the x contributions have aheady been taken 
into account in the previous section. 

We have restored the superscript ^'^' on the corrections of 5^1, in order to highhght the order in fi 
that we are working. 
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in the previous sections and in Appendices El |B] and O To begin with, we need the 
defining equations for ti 2 and TJ^^ given by (12. lip and (13.141) respectively but with 
the generahzed 61^2 (see (15.51) ) instead (of (15.71) ) while identity (13. lip is very useful. 
We also need the value of the integral J' defined in (15. 8p and given by (15. 9p . f lC.4p 
and table IC.ll as well as (13.170 and (jA.ip that define ri^2 and r, t] respectively. The 
final formula for gj^J is eventually given below. 



The Formula for gj^J 



/lU 



^S = . lim I -^fi'0{x^)0{x-)\ log (k\h - 61I) a.+ ( 2^2 ^,2 ' 
61.2^ (±6,0) V2 V / \r{^2[x^Y 



yj - ri, 



+ 



1 



2x+ 



fell - h 

\h2-hi\ 
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^2 _ ^2 



^-^i)^^M^ V2^2(xV )"^^^^^^ 



5f4__ lim l-^^9{x^)9{x )sech?7< — log ( /i;|62 — ^i| ) ^(t — ti^ 
bi,2->(±6,o) I 2 I 2r 



+ 



Cb^i - ^ (^sech2r/ + ^r9. (^^9.^^ J(ri,r2,r) + (l f^ 2 



^Vi = ^ lim <^ -^y.^e{x+)e{x- 



h 



1 — (^21 



ri 



bi,2-^(±b,0) I a/2 



|52-6i|2r2 + 2(x 



±^2 



^(r-n) 



2(^^21; 



rrrf^T sech^r? 



J{ri,r2,T) 



(5.10a) 



(5.10b) 



(5.10c) 



(2) r 

g\^ = hm 

bl,2^(±fe,0) 



5^!? = ^ lim 



-n'^6{x^)6{x )sechr7 



X {C62. +C +C.2.}^(^i'^2,r)l (5.10d) 



-IJ?0{x^)6{x )sechr7 
bi,2->(±b,o) I 4 



{^Lfen +C.2. +25L..2 +2Cfe22}^(^i'^2,r)l. (5.10e) 
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In order to arrive to f IS.lOp we have convoluted (15.11) with the right hand side of 
gSD and employed (Q), (COj) . f EOj) . (1X3|) . (COj) and (IRSJ) . In particular we have 
applied ( lAlej) and dEl for ( IQij) . ( 1X2]) and (lAlSJ) for (lilObj) . (IXil) and ( IXSJ) for 
(irOcj) and ([Oj) for both KMt and fHlQej) . 

The reason in preferring to work with the generahzed 61^2 is because the g 2 
may be obtain from g ^ under 1 4-> 2 before taking the limits as in (15. 7p : thus 
reducing the amount of calculations. Finally, (— /u) are obtained from (+/i) under 
the (simultaneous) interchanges (+ o — ) and (6i o 62)- These steps complete the 
determination of g^^J. Formula (IS.lOp is the final result of this project and we will 
analyze it in the next Chapter. 
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CHAPTER 6 
AREA OF VALIDITY AND CONCLUSIONS 



6.1 Area of Validity 

We have seen below (12. 4p that the parameter /i we use in our expansion has dimensions 
of length. On the other hand we know that the components of the metric should be 
dimensionless. Hence, each power in /i is compensated by an inverse power of the 
coordinates times a logarithm (with argument fcri 2, kb or kr) in some power at most. 
For simplicity (although not necessary), we restrict our discussion at middle -rapidity 
where ?] f^ that is for x^ ~ x~ (see flA.ip and figure IX!T|) . This means that any 
fc^'^-order (in yu) contribution to the metric, where fc is a positive integer, will generally 
have the form 

gl'H^n = ^i'e{x+)e{x-)J2c^Jk{xn-^ ^ + J + l = k>o (e.i) 

iji i ^ 

(see (I5.10p )where Cijk{x^) are dimensionless real functions of x^ Q We thus believe 
that our expansion is valid at high energies, that is small yU, compared to ri^2 and r. 
In fact, only one of ri^2 or t has to be large compared to yU. So for instance (IS.lOp is 
a good approximation for small r but large ri.2 (see figure 16.11 region /) and also for 



""^This form applies away from the light-cone because on the light-cone there also exist Sir — ri^2) 
terms (see (|6.2p ). 

^Such as 0-functions, logarithms and real coefficients 
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/'. T < r2 < ri 



I. T < ri < ra 




i^ X 



Figure 6.1: The reaction plane: Regions /// and ///' correspond to ri < r2 < r and 
r2 < ri < r respectively. The dark dots are the centers of the Shockwaves 
and are located at an impact parameter 2b apart while ri^2 denote the 
distance of the arbitrary point r from the center of each shockwave (right 
and left respectively). According to causality, at any given proper time r, 
the propagation from the centers will reach the points on the peripheries 
(at most). This suggests that any given point r on the transverse plane 
of the "produced" metric at given r will evolve according to the region 
where it belongs to (see equations (IS.lOp . flC.4|) and flC.5|) ) and there are 
three different possibilities (and three from their mirror images). 



small ri^2 but large r (region ///) provided that the massless particles creating the 
Shockwaves are not energetically enough (/i is small). This is in contrast to [12] where 
/i has dimensions of length to the negative third power and hence the expansion there 
was valid for early proper times. 

6.2 Conclussions 



In this work we have found the first non-trivial causal corrections to the problem of 
Shockwaves colhsions in gravity created by boosting two black holes to the speed of 
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light. The colhsion is assumed asymmetric and occurs at low energies (see previous 
section). In terms of Feynman diagrams, our result, formula (15.101) . corresponds to 
the resummation of the diagrams of figures 12.31 and 13. 2[ 

Our concussions are summarized as follows. 

1. The corrections to g^^, evolve non trivially and constrained by causality in 
an intuitive way. In particular, the behavior of g^^, at any point on the transverse 
plane, is determined from whether the propagation from the center of each individual 
nucleus has enough proper time to reach the point under consideration or not. Figure 
16. H is a snapshot taken at given proper time r and depicts the six kinematical regions 
where g^i, evolves differently while it has the general form 

9^'J = Oir2 - T)e{r, - r)4,(x^ b) + ^(r - n)e{r - t^)A'JJ{x^, b) 
+ {e{r2 - T)e{T - rl)41(x^ 6) + (6 o -b) } 
+ terms proportional to S{r — ri^2)- (6-2) 

The indices J, //, /// on A^,^ correspond to the regions /, //, /// of figure 16.11 

P 

respectively q The 5-function terms arise from differentiating the ^-functions of the 

right hand sides of fIS.lOp (see also (15. 9p ). They represent two Shockwaves centered 
at the center of each shockwave and expanding on the transverse plane with speed 
r (that is with the speed of light). This particular behavior of the metric was our 
initial motivation for dealing with this problem and our calculations confirm our 
earlier conjecture [T2] . 



^The (6 o —6) terms cover region I' ,11' and ///' since under this interchange we have ri o r2- 
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2. The presence of matter, T^^, and the back- reactions affect the metric (see for 
example fIS.lOap and f IS.lOcp ) not only on the forward light-cone but also inside. This 
implies that we cannot in principle solve Einstein's equations in vacuum (ignoring 
the point-like particles that create the shocks) arguing that we are away from the 
sources unless we know the boundary conditions that these sources enforce on the 
metric inside the light-cone. This is in analogy to classical electrodynamics: solving 
Laplace equation for the scalar potential away from a point charge sitting at the 
origin without specifying the boundary conditions, one may obtain the trivial (zero) 
solution which obviously is not the correct one. 

3. The presence of the impact parameter 6 is a necessary requirement and not an 
additional complication introduced in the problem. Mathematically this is obvious 
from the fact that both T^,^ and g^,^ diverge when the impact parameter b tends zero. 
From equation 13.201 we have seen that conservation of T^^^ is violated violently and 
behaves as V^T^y ~ 5±v6{0) ^ 0. The metric tensor also exhibits a problematic 
behavior in the zero impact parameter limit: the the formula for (yf^^-*, equation (15. 4p . 
diverges logarithmically when 6 —t- as it is evident (for instance) from equation 
(IS.lOap . A similar ultraviolet (UV) divergence appears in perturbation theory [29] of 
gauge theories. This suggests that a head on collision may not be investigated using 
classical gravity. Instead, one has to apply a quantum theory of gravity in the same 
way one can not predict the electron-positron annihilation (in head on collisions) 
using Maxwell's equations. One has to turn into Quantum Electrodynamics in order 
to describe the process and predict the production of two photons. One may argue 
that a black hole may be formed and hence hide the violation of conservation behind 
the event horizon. 
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4. For future projects we propose that one could plot g^^} for several (fixed) impact 
parameters as a function of r, x^ and x^ (at central rapidities where x'^ = x~) in 
order to visualize the evolution of the metric. A very important aspect we ignored 
in our analysis is the role of the ultraviolet cutoff k (see (I2.4p ) which (for the case 
of a single shockwave) seems to define an ergoregionu (with radius r^, = 1/k) on the 
transverse plane. It would be interesting to check how our solution gets modified 
for impact parameters b < 1/k which would imply that the two ergoregions of the 
Shockwaves overlap. Finally, one could compute the gravitational radiation, take the 
limit of 6 — > and compare the result with the one obtained by [H El E]- 



We thank Samir Mathur for a related and informative discussion. 
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Appendix A 
INTEGRATION OVER THE LIGHT-CONE PLANE 

In this appendix we perform the x^ part of the integrations resulting when the Green's 
function fl5.ll) acts on the right hand side of fl4.9p . Proceeding as in fl5.3p we find out 
that we have to deal with five different cases (there also exists a sixth case that is 
calculated in appendix iBj) . 

Case I: 6{x^)6{x^) terms 

This case is trivial and almost all terms of (14. 9 p behave in this way. Defining 

1, ,x+. ^ 1 



T = V2x+x- r] = -log( — ) : x^ = ^re^"^ (A.l) 

where r is the proper time and t] the rapidity (see figure lA?T] for a geometrical mean- 
ing) we have that 

G®6{x'^)S{x'-)fi{r') = —9{x+)9{x-)sechr] [ dV'6 {r - r') fi{r + P) (A.2) 

Attt J 

where we have shifted the integration variable setting r' ^ r' + f. 

Case II: 6'{x~^)6'{x~) terms 

This case is more complicated as we have to integrate by parts the 5-functions. 
Three kind of terms will appear: (a) Terms that differentiate the 9{x' — x ) terms 
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'? r3>T2> Tt x^ 






oo 



Figure A.l: The kinematical variables r and rj. The hyperbolas indicate curves of 
constant r and increase along x°. The straight lines are lines of constant 
rj and increase from left to right as the arrows indicate. Along x^ we 
have (r = 0, r; = ±00) while along x^ we have 77 = and x"*" = x~ . 



of (15.11) and hence produce 5{x — x ) (terms). But the presence 5{x' — x' ) 
forces the (5-function term appearing in (15. ip to become 5{—\f— r'\) which is zero. 
Hence these terms do not contribute, (b) We have terms that differentiate the 
denominator and these contribute to the integrations, (c) Finally, we have terms 
that either differentiate the (5-function of (15. ip only or both, the ^-function and the 
denominator. In order to evaluate these terms we exchange the d^'^ that act on 
5 I \/2{x^ — x'^){x^ — x'^) — \r — r'\\ with — 9^±. Shifting the transverse variable 
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as in the previous case and performing the x'^ integrations we find that the contri- 
bution of both the (b) and (c) terms is 

G^6'{x'+)6'{x'-)fjj{r') = --^e{x+)e{x-)sechr] 



1,2 1 o /lo 

sech rj -\ — tOj. —Or 
2 \T 



T 



£r'5{T-r')fn{f+r') 
(A.3) 

where the differential operator ^-sech^?] H — rdr I —d^ \ acts on the integral while d^- 

T^ 2 \T J 

denotes a partial differentiation with respect to r. 

Case III: 5'{x^)5{x^) terms 

This is a simpler version of case // and working in a similar fashion yields 

1 1 



G®6'{x'^)6{x'^)fni{r') 



27VTy/2 



e{x+)9{x^ 



-dr sech^T] 



l + e±2'? 2r 



(fPS (r - r') fin{r + ?). 

(A.4) 



Case IV: 6{x )6{x^) terms 



G<^6{x'^)e{x'^)fyi{?) = -^e{x^) 

An 



6 f^2x±(xT -x'^) - r'^ 
X / dx'^e{x'^)e{x^ - x'^) I d^r'^-^y-^-^-^ -;^W^fiv{r + r') 



V2 



{x^ + (xT - x'^)) 



j(^2x^{x^ -x'^) -r' 



2n 



r 



l-V29ix+)9ix-) I d'r'9ir-r')-j^^—-^frv{r + r'). 



(A.5) 
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Case V: S'{x^)e{x^)S^'^\f - h^) terms 

This is a combination of cases /// and IV with fi = f — hi with i = 1,2 as in 
(13171) . We have 



G®5'{x'^)e{x'^)fv{r')5^^\r- h) = - l ^^^-^0(x± - x'^)9{x'^ - x'^)5{x'^)e{x'^] 

4:71 



oo 



^(5 ('v/2(x±^^x^±)(x^3x^ 






fvih)^, +^^ \ r ^ ,^ eix^ - x'^)eix'^) / ^— 






(A.6) 



where in the first equahty we ignored a term similar to case // (see term (a)) when 
integrating by parts the 5'{x' ) while in the second equality we performed the 5{x' ) 
integration and exchanged d^i± with — i9j.±. The rest two steps are obvious. 
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Appendix B 
EVALUATING THE INTEGRAL (5.4) 

We wish to evaluate the expression (15.41) by performing the integration on both the 
hght-cone and the transverse plane. We begin by performing the x~^ and x~ integra- 
tions. Usiner that / dx~ / dx~6(x^) = x^9(x^) we find 



= fi^9(x+) / dx'~x'-e(x'-)9(x- -x'-) 



j[^/2^^ 



— -^e{x+)e{x-) [ d^r'\f-r'\ ^ -\^-^'\ ^(^_|^_f/|) 

271 y/2x+ J If -r'P + 2(x+)2 



X 



(^2,x'iV_Ltl,i.'i + t2,x'2V_Ltl,x''2) (B.l) 



where after the first equality we assume that the ti,2 = log{k\f— 61,2!) while the x 
dependence is displayed explicitly and it is integrated out after the second equality. 
The next step is to perform the transverse integrations. The trick here is to integrate 
by parts the V^ti 3.1,2 = 9^.1,2 (V^ti) terms. The by parts integration produces two 
kind of terms: (a) those that do not act on ^2,2:1,2 and (b) those that act on ^2,1:1, 2. 
But the terms of case (b) are proportional to 

Vli2Viti ~ 6{f- b2)S{f- 61) ~ (5(62 - 61) (B.2) 
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which is zero for non zero impact parameter 62 — &i while for zero impact parameter 
it diverges violently; we conclude that an impact parameter is necessary (see section 
16. 2p . We now proceed to the remaining terms. Exchanging 5a;/i,2 with —8x^,2 and 
using (12. 8 P equation (jB.ip gives 



G^ldx-jd.-{H„Vlt,„+t,^Vlt,y. 



2 I /■ 2 I -* ~/|2 

V2x+ [J |f-r'|2 + 2(x+)2 

X C^^^S{f' - hi) + (1^2) 



\r' — 65 



12 






(B.3) 



where 61,2 are given by (15.51) . 
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Appendix C 
EVALUATING THE INTEGRAL (5.8) 

We wish to calculate the integral ( 15 .Sp that we encountered in section 15.31 We have 
jr(ri,r2,r) = - — / / r'dr'd(f)'6{T - r')\og{k\r' + fl\)\og{k\r' + f*2\)- (C.l) 

^'^T Jo Jo 

The quantities ri 2 are given by fl3.17p . The trick here is to expand the logarithms in 

their Fourier space: log(A;r) = — / — with k serving as an ultraviolet cutoff. 

Expanding both logarithms and performing the angular integration one obtains 

Performing the trivial radial integration we find thai J' is now defined by 

Next we perform the / and q integrations. These integrals have been calculated in 
P^ : we summarize the procedure: In order to perform these integrations one has to 
expand ^0(^1^+ /|) in an infinite sum of products of the form Jn{T\q\) Jn{T\l\) with n 
an integer and do the angular integrals (of q and 1) first. This factors out the radial 
integrations over / and q into two independent integrals. Then one has to perform 
these integrations and finally sum over n. The final result reads 
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Table C.l: We defined r i 2 = r — 61 2 



cases 


e> 


V> 


^>V> 


^<V< 


J". 


region (see figure |fi.l ) 


1 


ri 


r2 


rir2 


r' 


Ji 


I 


2 


ri 


T 


rri 


rr2 


J2 


II' 


3 


r 


^2 


rr2 


rri 


Js 


II 


4 


T 


T 


r' 


rir2 


J4 


III 



Jir, n, ra) = ln(^>A;) ln{r]yk) + - 



L^, { e^»^) + LzJ e-'-^-^^ 



^>(<) = Tnax{min){ri,T) V>{<) = max{min){r2,T), 

(rl).(ri) = cos(a)rir2. 



(C.4a) 

(C.4b) 
(C.4c) 



So here a is tlie angle between fi and r2 and Li2 is the dilogarithm function while J^ 
is real as it should. Equation f lC.4|) implies that JT depends from the ordering of ri, 
r2 and r. There are in principle six distinct ways to order them. However it turns 
out that the cases ri^2 > t and ri^2 < t are interdependent from the relative ordering 
of ri and r2. This degeneracy reduces the possible cases to four which we organize 
by introducing table IC.lI rl that in turn help us to write a (unified) formula for J' 

J{ri,r2,T) = 9{ri - r)6'(r2 - r) Ji(ri,r2, r) + 6'(r - r2)6'(ri - r)^2(ri,r2, r) 

+ e^r - n)e{r2 - T)Jsir,, r2, r) + e{r - n)e{r - n)Mr,, r2, r). (C.5) 



^Js for instance means J'(C>^> — Tr2,£,<r]^ — rVi) with J'(C>'7>,'?<'7<) given from (|C.4I 
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